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Abstract  of  Dissertation  Presented  to  the  Graduate  Council 
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THE  ELECTROKINETIC  DETERMINATION  OF  THE 
STABILITY  OF  LAMINAR  FLOWS 

By 

Addison  Guy  Hardee,  Jr. 

August,  1968 

Chairman:   Knox  T.  Millsaps,  Ph.D. 

Major  Department:   Aerospace  Engineering  : 

An  experimental  investigation  of  the  transition  from  laminar  to 
turbulent  flow  in  a  pipe  of  circular  cross  section  has  been  conducted 
by  examining  the  electrokinetic  phenomena  generated  by  the  interaction 
of  the  fluid  flow  with  the  electric  double  layer.   It  was  found  that  ■ 
the  critical  Reynolds  number  of  the  flow  depends  on  the  electrolyte 
concentration  of  the  fluid,  that  Reynolds  number  being  higher  for  lower 
concentrations. 

The  investigation  also  revealed  that  the  fluid  flow  alters  the 
equilibrium,  state  of  the  electric  double  layer.   In  addition,  the  .. 
transition  to  turbulent  flow  can  be  detected  by  a  change  in  slope  of  the 
curve  for  streaming  potential  vs.  pressure,  but  only  if  the  concentra- 
tion of  the  electrolyte  is  a  very  low  value. 
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INTRODUCTION 


bcoDe 

The  present  investigation  is  the  natural  outgrowth  of  previous 
work  by  Anderson*'""',  of  which  it  represents,  at  the  same  time,  both 
an  extension  and  clarification.  Using  the  electrokinetic  phenomena 
associated  with  the  electric  double  layer,  the  transition  from  laminar 
to  turbulent  flow  in  a  pipe  is  studied  by  examining  the  growth  of 
perturbations  in  the  flow  very  close  to  the  pipe  wall.   The  electrokinetic 
phenomenon  of  interest  here  is  the  following:   when  an  electrolyte, 
however  weak,  flows  through  a  pipe,  a  potential  difference  can  be 
measured  between  a  pair  of  electrodes  placed  at  the  ends  of  the  pipe; 
this  is  the  streaming  potential  and  is  a  result  of  the  electric  double 
layer  which  exists  at  the  solid-liquid  interface  at  the  wall  of  the  pipe. 
The  electric  double  layer  is  discussed  more  fully  on  page  7. 

The  purpose  of  Anderson's  investigation  was  primarily  that  of 
determining  precisely  the  minimum  critical  Reynolds  number  of  Poiseuille 
flow  in  a  pipe;  thus  the  fluctuations,  which  he  discovered,  of  the 
streaming  potential  were  used  to  detect  the  onset  of  turbulence  in  the 
flow. 

The  radial  location  of  perturbations  in  the  flow  very  close  to. 
the  v?all  is  possible,  for  the  extension  of  the  electric  double  layer 


into  the  fluid  is  of  the  order  of  Angstroms  and  this  extension  (the 
double  layer  thickness)  is  a  function  of  the  concentration  of  the    , 
electrolyte.  The  present  investigation  examines  the  nature  of  the 
streaming  potential  fluctuations  for  double  layer  thicknesses  ranging 
from  539  to  about  5,000  A.  ' 

The  Stability  of  Laminar  Flows 

The  theory  of  the  stability  of  laminar  flows  is  usually  traced  to 

(2) 

Osborne  Reynolds    ,  who,  from  experimental  observation  and  theoretical  . 

studies,  postulated  that  the  state  of  laminar  flow  is  disrupted  by  the 
amplification  of  small  disturbances;  Reynolds,  however,  gives  credit  to 
Stokes  for  the  concept. 

Reynolds  treated  both  inviscid  and  viscous  fluid  stability;  many 
workers  since  have  spent  a  great  deal  of  effort  on  the  mathematics 
underlying  Reynolds'  hypothesis — namely,  to  superimpose  a  small  periodic 
perturbation  on  the  mean  flow  and  examine  the  growth  or  decay  of  this 
disturbance.   Most  notable  among  the  names  of  these  workers  are  those 

of  Rayleigh  "'    ,  Sommerfeld^   ,  Orr    ,  Tollmien    and  Prandtl 

(8)  (9)  ■ 

Mention  should  also  be  made  of  Heisenberg    and  Lorentz 

Just  one  example  of  Rayleigh 's  work  will  be  mentioned  here;  he 
showed  tha"C ,  for  frictionless  flow,  a  point  of  inflection  in  the  velocity 
profile  is  a  necessary  condition  for  instability,  i.e.,  for  amplifica- 
tion of  disturbances;  later  Tollmien  was  able  to  prove  that  this  is  also 
a  sufficient  condition. 

The  great  majority  of  theoretical  work  has  been  directed  toward 
the  stability  of  two-dimensional  flow,  the  application  of  the  theory 


being  the  superposition  of  a  two-dimensional  perturbation  on  the  stream 
function  of  the  mean  flow.   Neglecting  terms  quadratic  or  higher  in  the 
disturbance  leads  to  the  Sommerfeld-Orr  equation,  which  is  the  fundamental 
differential  equation  for  the  disturbance  amplitude. 

Use  of  a  two-dimensional  disturbance  was  questionable  until 
Squire ^■'"^^  showed  that  a  three-dimensional  disturbance  is  equivalent  to 
a  two-dimensional  one  at  a  lower  Reynolds  number — at  least  when  applied 
to  a  two-dimensional  flow;  this  has  become  known  as  Squire's  theorem. 

Heisenberg  was  the  first  to  deduce  theoretically  the  instability 
of  plane  Poiseuille  flow  for  sufficiently  large  Reynolds  numbers, 
although  he  did  not  calcu:.ate  a  lower  critical  Reynolds  number.   The 
theory  of  Heisenberg  was  extended  by  Tollmien     and  Schlichting    ; 

both  calculated  the  neutral  stability  curves  for  the  boundary  layer  over 

(13)  ■       • 

a  flat  plate.   Much  later,  Lin     performed  the  calculations  agaxn 

obtaining  more  clarity,  and  succeeded  in  calculating  the  curve  of  neutral 
disturbances  for  plane  Poiseuille  flow.   Contrary  to  the  conclusion  of 
Lin,  Pekeris     applied  a  different  technique  and  concluded  tnat 

Poiseuille  flow  is  stable  for  all  Reynolds  numbers.   To  resolve  this 

Cl5) 

disagreement,  Thomas     calculated  the  critical  Reynolds  number  by 

direct  numerical  methods  and  found  it  to  be  5780. 

The  more  difficult  analytical  problem  of  the  stability  of  pipe 
flow  has  not  been  resolved  with  such  clarity.  Sexl  was  the  first 
to  solve  the  viscous  problem,  although  only  for  axisymmetric  disturbances. 

The  results  have  been  questioned  because  of  his  mathematical  simplifica- 

(17) 

tions.   For  a  small  region  near  the  wall,  Pretsch     showed  the  problem 

became  the  same  as  that  of  a  disturbance  applied  to  plane  Couette  flow 


and  Pekeris     obtained  a  Eolutlon  for  the  region  near  the  axis  of  the 

( 18 ) 
pipe.   Corcos  and  Sellars     s^^^  s  solution  which  accounts  for  the   .  v 

work  of  both  Pretsch  and  Pekeris.   The  conclusion  dravjn  from  the  work  of 

these  investigators  is  that  Hagen-Poiseuille  flow  is  stable  for  small 

disturbances  and  the  work  of  Sexl  and  Spielberg     confirms  this.   (Sexl 

and  Spielberg  also  showed  that  Squire's  theorem  does  not  hold  for  axially 

symmetric  flows.)  Experimentally,  Leite     failed  to  observe  any 

amplification  of  small  axisymraetrical  disturbances  (placed  in  the  inlet, 

close  to  the  wall)  dov/nstream  in  a  circular  pipe  at  Reynolds  numbers  as 

high  as  13,000.  .   . 

All  of  the  investigations  to  date  imply  that  Poiseuille  flow  is 

stable  for  small  disturbances.   It  is  an  experimental  observation  that 

turbulent  flow  occurs  in  pipes — a  seeming  paradox  when  one  regards  the 

mathematical  results.   The  resclution  of  the  "paradox"  may  lie  in  the 

contrast  between  "small"  ar.a  finite  disturbances  occurring  in  the  flow 

or  in  the  symmetry  of  the  disturbance.   In  this  connection,  Meksyn  and 

(21) 
Stuart     showed  that,  in  a  channel  the  lower  critical  Reynolds  number 

decreased  as  the  amplitude  of  the  superimposed  oscillations  increased, 

which  is  in  accord  with  the  qualitative  observations  of  Reynolds.   A 

(22) 

possible  explanation  has  been  given  by  Gill    ,  who  reviewed  the  above 

theoretical  papers  and  indicated  questionable  steps  in  their  procedures. 

(23) 

On  the  other  hand,  Tatsumi     has  predicted  theoretically  that 

the  flow  in  the  inlet  of  a  pipe  is  unstable  at  a  Reynolds  number  of 
9.7  X  10  .   Experimentally,  both  Ekman  (unpublished)  and  Taylor 

.  Li 

(unpublished)  obtained  laminar  flow  up  to  Reynolds  numbers  of  5  X  10 

and  3.2  X  iO  ,  respectively;  Ekman  made  use  of  Reynolds  original  apparatus, 


The  Experimental  Investigations  of  Transition 

Hagen^"^^  first  noted  the  transition  from  laminar  to  turbulent 
flow  while  determining  the  law  of  resistance  for  pipe  flow.   He  was 
aware  that  the  "critical  point"  depended  on  the  velocity,  viscosity  and 
the  pipe  radius.   The  breakdown  of  what  we  today  call  laminar  flow  was 
noted  by  the  pulsing  of  the  jet  from  the  pipe  and  also  by  the  addition 
of  sawdust  to  the  flow,  showing  irregularities  present  above  the  critical 

point. 

The  fundamental  investigation  of  the  phenomenon  of  transition  was 
performed  by  Reynolds,  who  showed  conclusively  that  there  exist  two 

possible  modes  of  fluid  flow — laminar  and  turbulent .   He  was  most 

(25) 

probably  not  aware  of  Hagen's  work,  which  predated  that  of  Poxseuille^ 

On  the  other  hand,  Reynolds  was  in  possession  of  the  Navier-Stokes  equa- 
tions and,  by  dimensional  reasoning,  was  able  to  determine  the  form  of 
the  parameter  governing  the  "critical  point."  The  parameter,  of  course, 
is  the  Reynolds  number  and,  being  the  "similarity"  parameter  for  viscous 
flow,  is  more  than  just  the  parameter  for  transition.   Though  Reynolds' 
paper  is  often  quoted,  one  passage  from  his  1883  paper  is  worth  noting, 
especially  with  regard  to  the  aforementioned  "paradox."  Concerning  the 
sudden  disruption  of  the  flow,  he  writes: 

The  fact  that  the  steady  motion  breaks  down  suddenly 
shows  that  the  fluid  is  in  a  state  of  instability  for 
disturbances  of  the  magnitude  which  cause  it  to  break 
down.   But  the  fact  that  in  some  conditions  it  will 
break  down  for  smaller  disturbances  shows  that  there  is 
a   certain  residual  stability  so  long  as  the  disturbances 
do  not  exceed  a  given  amount . 

In  the  second  of  two  independent  experiments ,  Reynolds  determined  the 


minimum  critical  Reynolds  number  of  a  long  straight  pipe ;  the  value  he 
found  was  approximately  2000. 

Since  that  time,  many  workers  have  repeated  Reynolds'  experiments, 
some  with  interesting  variations.   The  most  extensive  repetition  was 
conducted  by  Stanton  and  Pannell^   •'.   Barnes  and  Coker^   '  used  a 
therm.al  m.ethod  of  detection  in  which  the  walls  of  the  pipe  were  heated 
and  the  onset  of  turbulence  was  detected  by  a  sharp  rise  in  the  tempera- 
ture of  the  interior  of  the  flow.   Reiss  and  Hanratty^^^^  developed  a 
technique  from  which  they  could  infer  the  behavior  within  the  so-called 
viscous  sublayer  by  measuring  the  mass  transfer  to  a  small  sink  at  the 
wall.   The  sink  was  a  polarized  electrode,  current  limited  by  mass 
transfer.  One  of  the  more  interesting  methods  is  that,  of  Lindgren    , 
who  made  the  disturbances  visible  by  using  polarized"  light  and  a  bi- 
refringent,  weak  solution  of  bentonite.   A  technique  utilizing  the 
electric  double  layer  was  developed  by  Anderson    and  used  to  determine 
the  lower  critical  Reynolds  number  for  Poiseuille  flow.   That  Reynolds 
number  was  found  ro  be  1907  ±3,  indicating  the  sensitivity  of  the 
technique. 

The  experiments  above  were  mainly  conducted  to  study  phenomena 
associated  with  "large"  disturbances,  while  it  was  not  until  Dryden's^'^^^ 
very  low  turbulence  wind  tunnel  became  available  that  experimentalists 

were  able  to  examine  the  small  disturbance  problem.   This  was  undertaken 

(31 ) 

by  Schubauer  and  Skramstad^   '^  who  made  an  almost  direct  transfer  of  the 

theoretical  method  to  the  physical  situation — oscillations  were  induced 
in  a  metal  ribbon  above  a  flat  place  placed  in  the  low  turbulence  tunnel 
and  the  amplification  of  the  flow  disturbances  measured  downstream  by 


means  of  a  hot  wire.   Their-  experiments  are  -.     regarded  as  excellent 
verification  of  stability  theory. 

The  experiment  of  Leite  '' '   mentioned  e.    .or  was  similar  to  that 
of  Schubauer  and  Skramstad,  but  axi-symmetric     urbances  were  super- 
imposed to  the  flow  of  air  in  the  inlet  of  a  pi^,  - 


The  Electric  Double  Layer 

The  phenomenon  known  as  the  electric  double  layer  has  been 
studied  extensively  by  chemists,  especially  in  connection  with  colloids 
and  with  electrode  processes. 

The  electric  double  lay^r  consists  of  an  excess  of  charge  present 
at  the  interface  between  two  phases,  such  as  a  solid  and  a  liquid,  and 
an  equivalent  amount  of  ionic  charge  of  opposite  sign  distributed  in  the 
solution  phase  near  the  interface.   Consider  one  phase  to  be  a  solid 
such  as  the  wall  of  a  pipe  and  the  other  to  be  a  weak  electrolytic 
solution. 

If  the  solution  is  caused  to  flow  past  the  wall,  such  as  in 
Poiseuille  flow  in  a  pipe,  there  develops  a  potential  difference  between 
the  ends  of  the  pipe  due  to  the  motion  of  the  distributed  charges.   This 
phenomenon,  known  as  the  streaming  potential,  was  discovered  by  Zollner^   \ 
and  subsequently  Helmholtz  "^  gave  an  explanation  based  on  Poiseuille 
flow  and  the  concept  of  the  double  layer  developed  by  Quincke    . 

The  early  workers  in  the  field  considered  the  double  layer  to  be 
composed  of  two  distinct  layers  of  charge,  one  fixed  to  the  wall  and  one 
free  to  move  with  the  fluid.   The  more  realistic  model  was  proposed  by 
Gouy  "■   and  Chapm.an     who  independently  formulated  the  theory  of  the 


diffuse  double  layer,  which  is,  in  essence,  the  theory  of  ionic  atmos- 

( 37 )        ( 38 ) 
pheres  given  some  ten  years  later  by  Debye  and  Huckel    .   Stern 

modified  this  theory  to  account  for  the  finite  size  of  the  ions  at  the 

wall. 

An  excellent  summary  of  the  classical  physics  of  the  effect  has 
been  given  by  Smoluchowski^"   and  extensive  analyses  of  the  approxima- 
tions used  in  the  various  theories  are  given  by  Kirkwood     and 
Casimir^^-^^ 

The  analysis  given  follows  Kruyt^    .   The  charge  at  the  interface 
is  considered  to  be  adsorbed  on  the  solid  surface  and  uniformly. distri- 
buted, while  the  solvent  is  assumed  to  be  a  continuous  media,  influencing 
the  double  layer  only  through  its  dielectric  constant.   Coulomb  inter- 
actions in  the  system  are  described  by  Poisson's  equation 

V  M'  ^  -  -^  (1)  . 

where  ^  is  the  potential  (having  a  value  of  lo   at  the  wall), 
^  is  the  charge  density, 
€  is  the  dielectric  constant,  ,  ■ 
and        V  is  the  Laplacian  operator. 
The  number  density  of  the  i   ionic  species  is  assumed  to  be  given  by 

where  ^io   is  the  number  density  of  the  i*-"  species  at  the  wall, 
2;  is  the  valence, 
C  is  the  electronic  charge, 
^  is  Toltzmann's  constant, 
and        "T  is  the  temperature  of  the  solution. 
The  charge  density  is  given  by 


Y   =  5.^-.e.n;  (3) 


Combining  equations  (1),  (2),  and  (3),  one  obtains 

i 

which  is  the  differential  equation  for  the  potential  as  a  function  of 
the  space  coordinates. 

Upon  assuming  an  infinite  plane  wall,  equation  (4)  is  simplified 
to 


I 


with  the  boundary  conditions 

where  u      is  the  distance  from  the  wall. 

The  first  integration  of  (5)  is  carried  out  after  multiplying 
through  by  2  <jwi  . 

Upon  integration  and  applying  the  boundary  conditions,  the  result 


is  obtained 


@'-  '-^In,.exF(-2,et/KT) 

The  equation  can  be  simplified  by  considering  a  single  binary 
electrolyte,  therefore 

(7) 

The  second  integration  is  performed  after  writing  equation  (7)  as 


10 


ScH'/K-i 


^^ __       \\[lM^ciM 


Therefore 


(8) 


where 

For  snail  potentials,  this  ir.ay  be  simplified  by  expanding  the  exponential 
terms  to  yield 


% 
5    -     —     f 


'X..  -   iy\ 


or 


H'  ^  %  eypC-'X^)  (9) 


showing  that  the  potential  decreases  exponentially  to  zero  over  a  distance 
of  the  order  of  magnitude  of  ^/y.    .      Thus  the  thickness  of  the  double 
layer  is  of  the  order  of  '/^  .   At  room  temperature  ^  is  approximately 
3  X  10  Z yfc'  where  C  is  the  concentration  in  gram -moles /liter.   As  an 
example,  the  double  layer  thickness  is  approximately  10"   cm  for  a 
O.OOIM  KCl  solution. 

The  Streaming  Potential 

The  streaming  potential  is  generated  when  the  electrolytic  fluid 
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is  caused  to  flow  by  applying  a  pressure  difference  between  the  ends  of 
the  pipe;  the  flow  displaces  the  charges  in  the  movable  portion  of  the 
double  layer  and,  so,  constitutes  an  electric  convection  current.  As  a 
consequence  of  this  current,  a  potential  difference,  the  streaming  '    '  ■ 
potential,  arises  between  the  ends  of  the  pipe.   In  the  equilibrium  state, 
an  equal  and  opposite  "conduction"  current  counterbalances  the  convection 
current;  the  conductance  determining  the  conduction  current  is  usually 
assumed  to  be  the  bulk  conductivity  of  the  fluid.  •    . 

Assuming  the  charge  density  on  the  wall  is  independent  of  the 
flow,  the  convection  current  is  ,.  ,- 

I.  =  ^^^^^  (10) 

where  U  is  the  hydrodynamic  velocity  and  A  is  the  cross-sectional  ' 


area. 


Substituting  from  Poisson's  equation  yields: 


Assuming  that  the  double  layer  extends  a  distance  out  from  the  wall 
which  is  small  compared  to  the  radius  of  the  pipe,  equation  (11)  may  be 
written: 

(12) 

Successive  integration  by  parts  and  application  of  the  boundary  condi- 
tions 


'•^ 
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and 


yields 


\\  =  o         «."t  ^   -    O 


^ 


-^\^i^);f-m'''^^h' 


(13) 


where  S  is  the  circumference  of  the  pipe  $nd  %  is  the  value  of  ^ 
at  y=0. 

If  the  flow  is  laminar,  '  _: ' 

/c)w.\     _  Ij:^ 


^<i:)i^.o    ^ 


where  '*^^^J  is  the  shearing  stress  at  the  wall  and  ^/o>-  is  the  viscosity 

of  the  fluid.   Also,  if  the  flow  is  in  a  constant  area  pipe, 

if 


where  X  is  the  coordinate  along  the  axis  of  the  pipe.   Therefore, 
equation  (11)  becones 


The  second  term  in  equation  (14)  will  be  neglected  since  it  may  be  shown 
that  the  ratio  of  the  second  rerm  to  the  first  is  the  order  of 
where  r  is  the  radius  of  the  pipe.   Therefore,  if  the  pressure  gradient 
is  constant  over  the  cross-section,  equation  (14)  may  be  written  as 

'i-^T^dx^  (15) 
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The  conduction  current  is  given  by  . 

ll--\\'^dX  (16) 


where  %     is  the  specific  conductance  of  the  solution  and  E.   is  the 
streaming  potential. 

If  "X  and  -j^        are  constant  over  the  cross-section,  equation 


(16)  becomes 


li  -    '^  dx  (17) 


At  equilibrium,  i,  =  -i^;  hence 


(18) 


Under  the  above  assumptions,  equation  (18)  shows  that  the  streaming 
potential  is  a  linear  function  of  the  pressure  drop  between  the  ends  of 
the  pipe. 

The  first  term  on  the  right-hand  side  of  equation  (13)  indicates 
that  the  convection  current,  and  hence  the  streaming  potential,  is  a 
function  of  the  velocity  gradient  at  the  wall. 


CHAPTER  II 


THE  EXPERIMENT 


General  Description 

The  present  investigation  was  undertaken  to  examine  more  closely 
several  aspects  of  streaming  potential  phenomena  which  the  experiments 
of  Anderson^  -^  brought  to  light. 

The  experiment  is  an  extension  of  Anderson's  work  in  which  he 
developed  a  technique  for  the  detection  of  turbulence;  therefore,  a 
short  account  of  his  experiment  and  the  reasoning  behind  it  will  be 
given. 

The  theory  of  Chapter  I  shows  that  the  streaming  potential  is 
directly  proportional  to  the  pressure  gradient  along  the  pipe.   Also, 
equation  (13)  shows  that  the  streaming  potential  varies  directly  and 
linearly  with  the  velocity  gradient  at  the  wall. 

If  we  now  restrict  the  discussion  to  Poiseuille  flow,  the  velocity 
distribution  across  the  pipe  is  given  by  ••  .  '  ; 


at  a  sufficient  distance  from  the  inlet,  where 

a  is  the  radius  of  the  pipe, 

?  is  the  pressure , 

X  is  the  coordinate  along  .the  pipe, 

14 


r-.   -.  -  .-/■: 
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__x/t  is  the  viscosity  of  the  fluid, 
r  is  the  distance  from  the  center  of  the  pipe, 
and        u  is  the  velocity. 

If  this  relation  is  substituted  into  equation  (13),  the  streaming 
potential  is  again  found  to  be  a  linear  function  of  the  pressure  gradient 
along  the  pipe;  this  has  been  observed  by  physical  chemists     for  many 
years.   Since  equation  (13)  indicates  the  dependence  of  the  streaming 
potenxial  on  the  velocity  gradient  at  the  wall,  Anderson  reasoned  that 
streaming  potential  measurements  could  be  used  to  determine  the  transi- 
tion from  laminar  to  turbulent  flow,  for  the  velocity  gradient  near  the 
wall  is  different  in  these  two  modes  of  flow. 

Anderson  looked  for  a  change  in  the  value  of  the  quantity  AE/AP 
(i.e.,  a  "break"  in  the  cur^ve  of  E  vs  P)  when  transition  occurred.   His 
experiment  was  conducted  using  an  electrolytic  solution  of  O.OOIN  KCl 
flowing  in  a  Pyrex  capillary  tube  having  a  diameter  of  0.0242  inches 
and  a  length  of  4.8  inches.   The  inlet  to  the  pipe  was  artificially 
roughened . 

Although  he  found  a  suggestion  of  such  a  break,  Anderson's  data 
did  not  unequivocably  show  one;  instead  he  discovered  unique  fluctua- 
tions in  the  streaming  potential  which  appeared  at  the  transition 
Reynolds  number.   Using  these  fluctuations  as  a  guide,  he  was  able  to 
determine  the  minimum  critical  Reynolds  number  very  precisely.   In 
addition,  he  showed  that  the  streaming  potential  did  not  vary  linearly  . 
with  pressure  gradient  if  the  stream.ing  potential  was  measured  across 
the  full  length  of  the  pipe.   This  is  the  effect  of  the  entrance  length 
(that  length  of  the  pipe  necessary  for  establishing  Poiseuille  flow) 
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which  the  theory  does  not  include .  Most  of  the  measurements  of  streaming 
potential  by  physical  chemists  have  not  taken  this  into  account. 

The  fluctuations  appeared  as  "spikes,"  which  always  represented 
an  increase  in  (positive)  voltage  near  the  transition;  above  transition 
the  spikes  appeared  to  be  positive  and  negative. 

The  present  investigation  was  undertaken  to  investigate  the  effect 
of  the  electric  double  layer  thickness  on  streaming  potential  measure- 
ments, the  thickness  being  a  function  of  the  electrolytic  concentration. 
The  electric  double  layer,  for  most  liquids  (including  the  electrolyte 
used  by  Anderson),  is  well  inside  the  so-called  "viscous"  sublayer; 
thus  the  mechanism  for  producing  the  streaming  potential  fluctuations 
was  not  apparent. 

Little  is  understood  of  the  actual  conditions  at  a  solid  surface 
in  fluid  flows;  in  a  liquid,  there  is  the  electric  double  layer  to  be 
considered.   The  electric  forces  in  the  double  layer  have  been  neglected 
in  solving  the  Navier-Stokes  equations  for  the  velocity  in  the  pipe 
(unless  they  are  contained  in  the  boundary  conditions).   In  addition, 
the  theory  given  for  the  electric  double  layer,  including  its  thickness, 
is  based  on  a  stationary  fluid — the  analysis  of  the  interaction  of  a 
flow  and  the  double  layer  is  complicated.   The  double  layer  extends  out 
from  the  wall  a  distance  on  the  order  of  Angstrom  units.  An  aqueous 
KCl  solution  of  O.OOIN  gives  rise  to  a  double  layer  thickness  of  95. 9A; 
this,  of  course,  is  the  distance  from  the  wall  to  the  "center  of  gravity" 
of  the  double  layer— the  double  layer  can  extend  out  much  further.   Such 
small  distances  from  the  wall  show  that  the  double  layer  is  within  the 
viscous  sublayer.   The  measurements  of  Reichardt     and  Laufer^'^^^ 
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seem  to  point  conclusively  to  the  existence  of  such  a  region  near  the  wall 
in  which  viscous  stresses  greatly  outweigh  inertial  stresses.   In  this 
light,  the  explanation  of  the  streaming  potential  fluctuations  becomes 
more  difficult.   Reichardt     predicted  the  fluctuations,  but  implied 
that  they  were  caused  by  perturbations  in  the  mean  pressure  gradient 
along  the  pipe — exactly  what  this  means  at  a  specific  point  very  near 
the  wall  is  unclear.   Reichardt  searched  for  the  fluctuations,  but  was 
unable  to  detect  them,  using  a  quadrant  electromieter .   Also  he  found  that 
the  value  of  AG/AP  was  essentially  the  same  for  both  laminar  and  turbulent 
flow. 

It  was  decided  to  extend  Anderson's  investigation  using  a  longer 
pipe  and  a  higher  gain  electrometer-amplifier  and  placing  greater 
emphasis  on  noise  reduction.   In  addition,  five  different  double  layer 
thicknesses  (ranging  from  589A  to  about  5,000A)  would  be  used,,  so  that 
the  double  layer  would  extend  over  varying  distances  into  the  viscous 
sub-layer.   The  longer  pipe  (543  diameters)  was  used  to  reduce  the  effect 
of  the  entrance  length. 

Laboratory 

The  laboratory  was  especially  designed  and  constructed  by  Anderson 
for  streaming  potential  experiments.   This  was  necessary  since  the  work 
is  particularly  sensitive  to  temperature  and  humidity.   The  thermal 
dependency  of  the  viscosity  of  the  fluid  necessitated  the  fine  control 
of  the  temperature,  v;hile  the  humidity  of  the  room  atmosphere  was  held 

very  low  due  to  the  nature  of  the  electrical  m.easurements .   The  currents 

_p      -2.0 
mvoive.i  were  of  the  order  of  10   to  10    amperes  with  a  source 
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9      11 
resistance  of  the  order  of  10  to  10   ohms;  a  high  humidity  level  could 

cause  surface  leakage  due  to  adsorbed  water  on  the  exterior  surfaces  of 

the  glass  apparatus . 

The  laboratory  is  a  10  by  16  foot  room  inside  an  air-conditioned 
building  and  the  entire  room  is  vapor  sealed  with  10  mil  polyethylene. 
The  floors  and  walls  are  insulated  with  3  inches  and  the  ceiling  with  4 
inches  of  Styrofoam.   Temperature  control  is  achieved  by  a  special  air- 
conditioning  system  which  holds  the  room  temperature  changes  to  ±0.75°  F. 
In  addition,  the  experimental  flow  apparatus  was  enclosed  within  a  small 
Styrofoam  box.  With  this  addition,  temperature  drifts  are  less  than 
0.02°  C  per  hour. 

The  air  drying  is  accomplished  by  use  of  a  compression-cooling- 
expansion  cycle.   This  enables  the  room  air  to  be  held  below  20%  R.H. 

Apparatus 

The  apparatus  used  is  similar  to  that  used  by  Jones  and  Wood    , 

(47) 
Kruyt    ,  Anderson  and  others  and  is  shown  schematically  in  Fig.  1. 

The  reservoirs  R^  and  R^  are  5  liter  boiling  flasks  fitted  with  ground 

glass  taper  joints,  J,  and  J„.   Tubes  L,  and  L_  are  6mm  Pyrex  with 

taper  joints  at  the  top  for  inserting  electrodes.   A,  and  A„  are  1/4 

inch  tubes.   The  tube  T  is  a  precision  bore  capillary  tube  0.0531  ± 

.0001  inches  in  diameter  and  34.3  inches  long.   The  capillary  tube  is 

connected  to  A^^  and  A  by  flanges  incorporating  Teflon  gaskets.   In 

operation,  the  liquid  was  forced  up  L,  and  through  the  capillary  by 

applying  gas  pressure  through  V-,. 

The  electrodes  E^^  and  E2  are  Ag-AgCl  formed  on  spirals  of  24 
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AWG  platinum.   The  Ag-AgCl  coating  was  made  with  a  mixture  of  90%  silver 
oxide  and  10%  silver  chlorate  formed  into  a  paste,  applied  to  the  wire 
and  heated  at  500°C  for  fifteen  minutes.      ..•-■■. 

The  gas  pressure  was  supplied  by  commercially  bottled  nitrogen. 
The  pressure  was  measured  by  a  simple  mercury  manometer  used  in  combina- 
tion with  a  micro-telescope  mounted  on  a  vernier  height  gage.  The 
pressure  could  easily  be  controlled  to  within  ±0.005  inches  Hg.   The  two 
three-way  stopcocks  in  the  supply  line  provided  a  means  of  applying 
pressure  to  the  flasks  or  venting  them  to  the  atmosphere. 

Cleaning  of  the  Apparatus 

The  glassware  in  the  apparatus  was  soaked  in  chromic  acid  for 
several  hours,  cleaned  with  hot  chromic  acid  and  rinsed  in  hot  conductiv- 
ity water.   The  glassware  was  then  leached  in  conductivity  water  over- 
night.  The  electrodes  were  leached  in  conductivity  water  for  several 
days . 

The  conductivity  water  used  as  the  solution  was  twice-distilled, 
having  nitrogen  bubbled  through  it  after  each  distillation. 

Electrolytic  Solutions  Used  in  the  Experiment 

The  experiments  were  conducted  with  5  concentrations  of  KCl 
solutions,  beginning  with  conductivity  water  as  the  first  solution.   Salt 
was  added  to  this  to  form  the  subsequent  concentrations.   Table  1  gives 
the  values  of  the  five  concentrations  and  the  corresponding  thickness  of 
the  double  layer. 


TABLE  1 
THICKNESS  OF  THE  ELECTRIC  DOUBLE  LAYER 


Concentration 


Norri.ality 


6(cm)"' 


1 
2 
3 
4 
5 


"6  =  double  layer  thickness. 


The  Electrical  Measurements 


oo 


6(A)-- 


oo 


0.5354X10 

4.16X10"^ 

4,160 

0.8040X10"^ 

3.40X10"^ 

3,400 

1.877X10"^ 

2.23X10"^ 

2,230 

25.800X10'^ 

0.589X10"^ 

589 
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Three  types  of  electrical  measurements  were  made  for  each  concen- 
tration of  electrolyte.   The  magnitude  of  the  streaming  potential  and 
the  root-mean-square  value  of  the  streaming  potential  fluctuations  were 
measured  as  functions  of  pressure.   In  addition,  the  resistance  of  the 
electrolyte-pipe-electrode  "source"  was  measured  before  and  after  each 
day's  run. 

The  streaming  potential  magnitude  was  determined  by  the  system 
shown  in  Fig.  2.   Using  a  General  Radio  Model  1230-A  Electrometer  (input 
impedance  10^    )  as  a  null  indicator,  the  streaming  potential  was 
balanced  by  the  output  of  a  0-90V  two-stage  voltage  divider.   This  was 
then  reduced  by  a  factor  of  1000  by  the  Leeds  and  Northrup  Volt  Box  and 
read  on  a  Leeds  and  Northrup  millivolt  potentiometer.   Streaming  poten- 
tials for  the  1.377X10"^N  and  26.80X10"S  solutions  were  read  directly 
from  the  General  Radio  electrometer. 
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The  R.M.S.  values  of  the  streaming  potential  fluctuations  were 
determined  using  the  system  shown  in  Fig.  3. 

The  streaining  potential  fluctuations  were  amplified  by  a 
differential  electrometer-amplifier  specially  constructed  for  electro- 
chemical measurements;  the  input  stages  are  Philbrick  Model  SP2A 

14 
operational  amplifiers.  The  input  impedance  of  this  amplifier  is  10 

ohms  and  the  gain  was  8.80  throughout  the  experiment  (gain  can  be  set 

at  8. SO,  88.0,  or  830).   The  output  of  the  voltage  divider  was  used  to 

cancel  the  D.C.  component  of  the  streaming  potential  (this  ranged  up  to 

80V).   The  output  of  the  amplifier  was  fed  into  a  strip-chart  recorder 

for  a  visual  record  of  the  fluctuations  and  into  the  Flew  Corp.  Random 

Signal  Voltmeter  for  their  R.M.S.  value. 

The  electrical  system  represented  by  the  electric  double  layer 

is  very  sensitive  to  stray  capacitance — it  was  found  that  an  excessive 

amount  of  "grassy"  noise  (which  almost  obliterated  the  smaller  signals) 

was  generated  by  vibrations  in  the  screen  panels  of  the  screen-room 

(Faraday  cage).   This  was  corrected  by  glueing  aluminum  foil  shielding 

onto  the  Styrofoam  thermostat  in  which  the  flow  system  was  placed.   This 

shielding,  along  with  the  usual  care  in  reducing  vibrations  in  all 

electrical  components,  gave  a  great  reduction  in  noise  level.   Typical 

signal  levels  were :  '  .  '  '       '. 

spike  magnitude  (peak)       O.IV 

spike  magnitude  (R.M.S.)     0-50mv 

streaming  potential  (D.C.)    0-85V. 
A  typical  spike  trace  is  shown  in  Figure  4. 

As  previous  experimenters  had  noted,  the  values  of  the  streaming 
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Figiire  U.   Tracing  of  Typical  Spikes  (Reynolds  Number  is  2572, 
Concentration  No.  2) 
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potentials  could  not  be  reproduced  accurately  from  day  to  day.  For 
each  separate  electrolytic  solution,  the  streaming  potential  showed  a 
long-time  drift  of  a  period  ranging  from  a  few  days  to  weeks,  reaching 
a  minimum,  and  then  increasing  slowly.  The  D.C.  component  of  the 
streaming  potential,  for  a  given  electrolyte  concentration,  was  taken 
over  a  two-hour  period  on  the  day  the  potential  reached  a  minimum. 
Reproducibility  was  also  improved  if  the  fluid  was-  pumped  at  high  flow 
rates  before  taking  data;  also,  higher  Reynolds  number  measurements  were 
alternated  with  those  of  lower  Reynolds  number.   If  this  alternation 
were  not  carried  out,  the  D.C.  values  for  the  low  flow  rates  (below 
Reynolds  numbers  of  about  900)  would  drift  upward.   This  indicates  that 
the  equilibrium  state  of  the  double  layer  is  changed  when  the  electrolyte 
is  moving . 

Determination  of  Reynolds  Number 

All  of  the  data  was  taken  as  a  function  of  the  pressure  difference 
across  the  pipe.   The  Reynolds  number  was  determined  by  mass  flow  rate 
measurements  using  the  last  (highest  concentration)  solution.  Reynolds 
numbers  were  taken  from  this  data,  shown  in  Figures  5  and  6.   The 
resistance  coefficient  is  defined  by    , 

where  Uvn  is  the  mean  velocity, 

■  .  -'   Lp  is  the  length  of  the  pipe 
and  .   ^    yo  is  the  fluid  density. 
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'igure  6.   Resistance  Coefficient  vs.  Reynolds  Number 
(Log-Log  Plot) 


CHAPTER  III 
RESULTS  AND  CONCLUSIONS 
R.M.S.  Component  of  Streaming  Potential 

The  R.M.S.  values  of  the  fluctuations  in  the  streaming  potential 
are  shown  in  Figures  7  to  11,  plotted  against  the  Reynolds  number;  these 
values  have  been  divided  by  the  gain  of  the  electrometer-amplifier.   The 
points  denoted  by  triangles  represent  laminar  flow,  i.e.,  no  spikes 
(fluctuations)  on  the  streaming  potential;  those  denoted  by  circles 
represent  the  occurrence  of  spikes.   Below  transition  (no  spikes)  and 
above  a  Reynolds  number  of  2560,  all  of  the  data  are  shown  on  the  graphs; 
in  between,  only  points  selected  to  indicate  the  mean  curve  are  given 
to  avoid  crowding.   The  average  deviation  from  the  curve  is  about  ±7  in 
Reynolds  number.   In  Figure  11,  all  the  available  data  are  shown  on  the 
graph . 

Each  figure  displays  the  same  general  trend;  a  peak,  followed  by 
a  dip,  and  a  second,  higher  peak  followed  by  a  rapid  decrease.   Such 
regular  variations  are  not  detectable  in  mass  flow  rate  measurements 
through  the  transition  range  of  Reynolds  numbers — only  an  irregular 
increase  (plot  of  resistance  coefficient  vs.  R  -  Fig.  6).   At  the 
highest  Reynolds  number  of  the  experiment,  transition  to  fully-developed 
turbulence  was  not  complete  (deduced  from  Fig. 6). 

The  frequency  of  occurrence  of  the  spikes  at  their  first 
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appearance  was  taken  as  0.002  per  second  (one  spike  in  an  eight  minute 
period)  and  at  R^   =  2755  was  approximately  2.1  per  second;  this  frequency 
was  determined  by  manually  counting  the  peaks  of  the  spikes  over  an 
eight  minute  period  of  the  chart  from  the  recorder.   Thus,  although  the 
frequency  response  of  the  electrometer-amplifier  and  electrodes  is  flat 
for  this  range,  the  response  of  the  system  comprised  of  double  layer- 
electrodes-amplifier  may  be  responsible  for  the  variations  in  the  Ej^j^g 

vs.  R  curves.   The  period  of  an  individual  spike  was  about  2  seconds. 

n 

The  existence  of  some  fluctuation  below  the  appearance  of  the 
first  spikes  is  a  measure  of  "free-stream"  turbulence  as  well  as  fluctua- 
tions from  the  inlet  of  the  pipe.   The  corresponding  R.M.S.  values  at  a 
Reynolds  number  of  1000  was  about  half  that  at  2U00. 

The  curves  for  each  double  layer  thickness  are  shifted  progressively 
to  lower  Reynolds  number  as  the  double  layer  thickness  decreases .  The 
first  spikes  appear  at  R  =  2585  for  the  first  concentration  and  at 
R„  =  2^:70  for  the  last  concentration.   The  peaks  and  other  salient 
features  of  the  curves  display  corresponding,  although  smaller,  shifts  in 
Reynolds  number.   The  magnitudes  of  the  shifts  are  too  large  to  be 
explained  by  changes  in  average  viscosity  of  the  fluid,  but  the  same 
cannot  be  said  of  the  viscosity  in  the  region  of  the  double  layer  as 
the  double  layer  readjusts  to  an  overall  concentration  change.   As  will 
be  shown  in  the  next  section,  the  equilibrium  state  of  the  electric 
double  layer  is  changed  by  the  flow  impressed  upon  it — whether  this  is 
a  change  in  potential  distribution  alone  or  with  a  simultaneous  change 
in  potential  at  the  wall  can  not  be  determined.   The  part  the  electric 
double  layer  itself  plays  in  the  stability  of  laminar  flow  is  unknown. 
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although  some  insight  into  its  role  might  be  gained  by  dimensional 

analysis. 

Consider  the  two-dimensional  steady  flow  of  a  viscous  fluid 
having  an  electric  double  layer  present  at  the  solid  boundaries.   The 
Navier-Stokes  equations  are 

where  ?.  is  the  electrostatic  force  on  a  fluid  element.   If  one  now 

1 

considers  two  flows  with  geometrically  similar  boundaries  and  let  L^ 

be  any  length  of  the  first  flow  field  and  L^   the  corresponding  length  in 

the  second,  then  L^   =  C-^L-^.      Similarly,      ,: 


(2) 


The  individual  terms  for  the  second  flow  field  in  relation  to  the 
corresponding  terms  of  the  first  flow  field  are 

Nil  C.7.  c>^\ 

_t_  ^  _  _L_  52.  ^ 

f*  ^X^    C^P,  C,  hi^,   '  (5) 

and  "^^^  ""  ^^"^  <^  '^""  ^^^ 
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In  order  that  the  equations  of  motion  for  the  two  flow  fields  may  be 
identical,  the  following  must  be  true  . 


The  equality  ~=-  -     -—ic"  ™^y  ^®  written 


(7) 


(8) 


which  means  the  Reynolds  numbers  of  the  two  flows  must  be  equal  in  order 
for  the  flows  to  be  geometrically  similar. 

The  relation  -^  •=.  ^    ~  means  that 


must  also  be  true .   This  may  be  written  as 


(9) 


FL? 

=  constant  (10) 


where  R  is  the  Reynolds  number  and  F  is  the  electrostatic  force  per  unit 
volume  on  a  fluid  element. 

The  electrostatic  force  per  unit  volume  on  the  fluid  element  is 
given  by 


4± 


}-^-Y-r^  (11) 
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where  q  is  the  electric  charge  contained  in  the  element  and  H'  is  the 
(electric)  potential.   For  flow  within  the  electric  double  layer 

where   4  is  the  charge  density.   Therefore,  the  parameter  in  equation 
(10)  may  be  written  ,    = 

or,  approximately  (if  the  double  layer  thickness  is  taken  as  the 
characteristic  length)  as 

Therefore ,  the  parameter 

1. 


/<fo  L  g.yp(-') 

must  be  a  constant  for  dynamically  similar  flows. 
,  The  charge  density  is  given  by 


(13) 


(14) 


where  n^^  is  the  number  density  of  the  negative  ions  and  x\^   is  the  number 
density  of  the  positive  ions.  Again  the  number  density  is  assumed  to 
be  given  by  a  Boltzmann  distribution: 

■   Ai=  n;,expC-Z;e.4'/KT^  .  (15) 

where  n^^  is  the  number  density  of  the  i   ionic  species  at  the  wall 
(the  wall  is  considered  to  have  a  negative  charge).'  Therefore  (for  a 
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single  valence,  binary  electrolyte),  the  charge  density  is  approximately 
given  by 


y.-e^.,e.p[-'^--^-P^-^^^] 


The  similarity  parameter  now  becomes 


-  e  4-,  L  no,  ^Y.Vl KT  1  «-^?^->^ 


^<^v  R, 


(16) 


The  quantities  4o  »  L  and  the  exponential  term  in  the  numerator  all 
increase  as  the  bulk  concentration  of  the  electrolyte  decreases;  there- 
fore ,  the  Reynolds  number  must  increase  in  order  that  the  parameter 
remain  constant.   For  more  dilute  solutions  (i.e.,  thicker  double  layer), 
the  higher  must  be  the  Reynolds  number  of  flow  phenomena  within  the 
electric  double  layer.   The  Reynolds  number  here,  of  course,  is  based 
on  the  thickness  of  the  electric  double  layer  and  the  local  velocity  of 
the  flow. 

The  relation  of  the  corresponding  mean  Reynolds  number  of  the  flow 
to  a  local  Reynolds  number  within  the  double  layer  depends  on  the  rela- 
tion of  their  respective  velocities;  whatever  this  exact  relation  may  be, 
the  parameter  given  by  the  above  analysis  agrees  qualitatively  with  the 
experimental  results — the  transition  to  turbulence,  as  indicated  by  the 
double  layer,  occurs  at  higher  Reynolds  numbers  for  more  dilute  solutions. 


Streaming  Potential 


The  values  of  the  D.C.  com.ponent  of  the  streaming  potential  for 
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the  five  solutions  are  shown  as  functions  of  pressure  (non-dimensional) 
in  Figures  12  to  17.  Excepting  a  region  at  the  beginning  of  transition, 
the  slope  for  laminar  flow  is  different  than  the  slope  for  turbulent 
flow.   Thus,  the  change  in  slope  can  be  used  to  detect  the  transition 
from  laminar  to  turbulent  flow.   A  comparison  of  the  graph  of  E^^  vs. 
pressure  for  the  last  concentration  (Fig.  17)  and  resistance  coefficient 
(?p)  vs.  Reynolds  number  (Fig.  6)  shows  that  the  streaming  potential 
indicaxes  a  slightly  higher  critical  Reynolds  number  than  mass  flow  rate 
measurements.   The  break  in  the  Fp  vs.  R^   curve  begins  when  the  eddies 
grow  enough  in  strength  and  number  to  increase  appreciably  the  resistance 
to  the  flow.   This  should  coincide  very  nearly  with  the  establishment 
of  a  different  velocity  profile  near  the  wall — thus  changing  the  value  of 

■  ^-  as  governed  by  equation  (13),  Chapter  I.  The  difference  in 
critical  Reynolds  number  given  by  the  two  methods,  mass  flow  rate  and 
streaming  potential,  may  indicate  that  the  turbulent  velocity  profile 
is  established  "close"  to  the  wall,  i.e.,  on  the  order  of  10  A,  at  a 
higher  Reynolds  number  than  shown  by  an  increase  in  flow  resistance.  A 
cuestion  arises  when  one  considers  that  the  velocity  gradient  of  equat 
(13)  m.ust  be  that  within  the  double  layer,  which  m.eans  well  within  the 
"viscous"  sublayer;  according  to  the  accepted  model  of  this  sublayer, 
turbulent  eddies  (necessary  to  produce  a  change  in  velocity  profile)  so 
near  the  wall  musx  be  very  weak.  Further,  it  is  usually  assumed  that 
the  velocity  profile  within  the  viscous  sublayer  does  not  differ  in 
laminar  and  turbulent  flow.  This  says  nothing  about  the  magnitude  of  a 
change  in  velocity  profile  which  the  double  layer  can  detect.  When 
'viscous"  sublayer,  it  should  be  kept  in  mind  that  the 
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Figure  13.   D.C.  Component  of  Streaming  Potential  vs.  Pressure, 
Concentration  No.  2. 
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Figure  14.   D.  C.  Component  of  Streaming  Potential  vs.  Pressure, 
Concentration  No.  3 
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entire  flow  in  the  pipe  is  "viscous."  Also  the  work  which  points  to  the 
existence  of  the  sublayer  is  based  largely  on  hot-wire  measurements  and 
these  have  only  been  made  over  the  outer  part  of  the  sublayer. 

Another  fact  bearing  on  the  above  question  is  that  in  deriving 
the  equations  for  the  streaming  potential,  it  was  assumed  that  the  flow 
conditions  (including  the  velocity  profile)  were  the  same  all  along  the 
length  of  the  pipe;  in  reality,  the  theoretical  conditions  are  attained 
only  after  the  entrance  length  of  the  pipe  has  been  traversed.   Further, 
the  entrance  length  necessary  for  establishment  of  such  conditions  is 
much  shorter  for  turbulent  than  for  laminar  flow.   This  means  that  the 
entrance  length  is  changing  in  the  transition  region  of  Reynolds  numbers 
(it  is  also  a  function  of  Reynolds  number  for  laminar  flow  and,  to  a 
lesser  degree,  for  turbulent  flow). 

The  velocity  profile  for  laminar  flow  in  the  entrance  length  is 
similar  to  that  for  established  turbulent  flow;  thus,  for  laminar  flow 
at  low  Reynolds  numbers ,  the  value  of  AE/aP  should  be  intermediate  to 
those  of  established  laminar  flow  and  established  turbulent  flow.  This 
is  seen  to  be  the  case  for  all  the  data  in  Figures  12  to  17 — the  slope 
of  the  laminar  portion  of  the  curve  decreases  with  increasing  Reynolds 
number. 

Measuring  the  streaming  potential  and  pressure  across  one  segment 
of  the  pipe,  excluding  the  entrance  length,  would  better  fulfill  the 
conditions  of  the  theory  and  show  whether  or  not  the  "non-developed" 
flow  in  the  changing  entrance  length  is  responsible  for  the  difference 
in  slopes.   This  was  considered  in  preliminary  planning  of  the  investiga- 
tion; however,  it  was  evident  that  any  method  of  isolating  one  length 
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of  the  pipe  would  introduce  disturbances  to  the  flow  and  would  require 
such  a  large  system  that  purity  of  the  solution  could  not  be  maintained 
with  certainty . 

The  magnitudes  of  streaming  potential  given  here  include  the 
average  value  of  the  streaming  potential  fluctuations.   These  average 
values  should  be  of  the  same  order  of  magnitude  as  the  R.M.S.  values, 
i.e.,  less  than  50  millivolts.   This  is  a  sufficiently  small  part  of  the 
Eq^  values  to  be  ignored.   The  data  for  the  last,  and  most  concentrated, 
solution  (Fig.  17)  shows  an  unexplained  peak  at  the  beginning  of  transi- 
tion.  The  peak  occurs  roughly  at  the  same  Reynolds  number  as  the  maximum 
R.M.S.  value  (Fig.  11),  but  is  about  twice  the  expected  magnitude. 
Figure  16  shows  Eq„  taken  from  a  Dymec  (Hewlett-Packard)  Integrating 
Digital  Voltmeter,  Model  2'401A,  set  up  as  shown  in  the  inset  to  this 
figure.   Integration  does  not  alter  the  basic  curve,  but  does  slightly 
change  the  ratio  of  the  turbulent  slope  to  the  laminar  slope.  These 
data  were  taken  with  concentration  4  on  the  same  day  as  the  data  shown 
in  Figure  10. 

Although  all  the  curves  are  similar,  they  represent  a  wide  range 
of  streaming  potential  magnitude.   In  Fig.  18  is  shown  how  the  ratio  of 
the  turbulent  slope  to  the  laminar  slope  varies  with  double  layer  thick- 
ness; the  slope  for  laminar  flow  was  taken  to  be  that  of  the  portion  of 
the  curve  immediately  preceding  transition.   This  ratio  approaches  unity 
for  a  sufficiently  small  double  layer  thickness. 

Conclusion 

The  conclusions  drawn  from  the  experimental  investigation  are 
delineated  as  follows. 
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^igure  18.   Ratio  of  Lair.inar  Slope  to  Turbulent  Slope  vs. 
Double  Layer  Thickness 
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The  critical  Reynolds  number  of  pipe  flow,  as  indicated  by- 
streaming  potential  fluctuations,  depends  on  the  thickness  of  the  electric 
double  layer  present  at  the  wall  of  the  pipe,  i.e.,  on  the  concentration 
of  the  electrolytric  fluid. 

The  equilibrium  state  of  the  electric  double  layer  in  a  stationary 
fluid  is  different  than  that  for  a  flowing  fluid,  at  least  for  laminar 
flow. 

The  transition  from  laminar  to  turbulent  flow  can  be  detected  by 
noting  a  change  in  slope  of  the  graph  of  E^q  vs.  pressure,  but  only  if 
the  electrolytic  solution  is  very  dilute. 
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